General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



COO-3077-149 


Courant Mathematics and 
Computing Laboratory 


U.S. Department of Energy 

(NASA-CR- 155925) A ViELL ECSFE BOON t A F Y N78-17343 

“ A L r l E PROBLEM IN 7PANJCN1C GAS DYNAHICo (Nev 
York Univ.) 48 p AC A03/HF AC 1 CSCL 20D 

One la a 

iiJ/34 C5690 


A Well Posed Boundary Value Problem 
in Transonic Gas Dynamics 

Jose M. Sanz 

l 


Research and Development Report 

Prepared under Contract EY-76-C-02-3077 
with the Office of Energy Research, NASA- Langley 
Research Center Grant NCR- 33-01 6- 167 and 
NASA-Ames Research Center Grant NGR-33-016-201 . 


Mathematics and Computing 
sbruary 1978 



New York University 



UNCLASS I F I ED 


Courant Mathematics and Computing Laboratory 
New York University 


Mathematics and Computing COO-3077-149 

A WELL POSED BOUNDARY VALUE PROBLEM IN TRANSONIC GAS DYNAMICS 

Jose M. Sanz 


February 1978 


Prepared under Contract EY-76-C-02-3077 with the Office of 
Energy Research, NASA-Langley Research Center Grant NGR-33- 
016-167 and NASA-Ames Research Center Grant NGR-33-016-201 . 


UNCLASSIFIED 



TABLE OF CONTENTS 


1. INTRODUCTION 1 

2. AN APPLICATION 5 

2.1 Equations of Motion and the Hodograph 

Transformation 5 

2.2 The Inverse Problem for Supercritical 

Airfoil Design 10 

3. GREEN'S AND NEUMANN'S FUNCTIONS AND THE 

POISSON FORMULAS 14 

3.1 The Fundamental Solution 15 

3.2 The Green's Function for the Unit Semicircle 18 

3.3 The Poisson Formula for the Unit Semicircle 21 

3.4 The Mixed Dir ichlet-Neumann Problem 25 

4. THE REFLECTION PRINCIPLE 30 

5. SOLUTION OF THE BOUNDARY VALUE PROBLEM 33 

BIBLIOGRAPHY 39 

FIGURE 43 


iii 


ABSTRACT 


A boundary value problem for the Tricomi equa- 
tion is studied in connection with transonic gas dynamics. 

The transformed equation 

Au ♦ yj- Uy * 0 

in canonical coordinates is considered in the complex domain 
of two independent complex variables. A boundary value 
problem is then set by prescribing the real part of the 
solution on the boundary of the real unit circle. 

The Dirichlet problem in the upper unit semicircle with 
vanishing values of the solution at Y ■ 0 is solved explicitly 
in terms of the hypergeometric function for the more general 
Euler-Poisson-Darboux equation. An explicit representation 
of the solution is also given for a mixed Dirichlet and 
Neumann problem for the same equation and domain. 

Reflection rules are given for these solutions which 
permit one to extend them from the upper to the lower unit 
semicircle. The transonic boundary value problem is solved 
by expressing the solution as a linear combination of these 
two types of solutions. 
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1 . 


INTRODUCTION 


The Tricomi equation 

(1.1) yu + u - 0 

1 xx yy 

is an example of a partial differential equation of mixed 
type, that is, an equation which is partly elliptic and 
partly hyperbolic. This equation has been widely studied 
because of its relation to transonic gas dynamics 110]. 

The Tricomi problem is a model of a well posed problem 
for equation (1.1). We recall [10] that it is formulated 
in a domain which is intersected by the parabolic line 
y = 0. In the elliptic half-plane the domain is bounded 
by a simple curve ending at two points on the parabolic 
line. In the hyperbolic half-plane the domain is bounded 
by two intersecting characteristics that start at these 
two points of the parabolic line. The Tricomi problem is 
well posed because the boundary values of the solution are 
prescribed on the whole elliptic boundary but only on 
one of these characteristics. 

An important problem in transonic gas dynamics is the 
computation of the flow past a given airfoil with speed at 
infinity close to but smaller than the speed of sound. In 
this situation the flow over the upper side of the airfoil 
may become supersonic. The partial differential equations 
of motion are then of mixed elliptic and hyperbolic type. 
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The problem of calculating continuous transonic flow past 
such an airfoil was shown by C. Morawetz ( 20 ] to be not 
well posed. 

Motivated by the method of complex characteristics, a 
transonic boundary value problem was developed in [9] for a 
certain domain obtained by analytic extension to complex 
values of the independent variables. Complex conjugate 
characteristic coordinates £ and n were introduced. In 
these coordinates the stream function iM£,£) is real in the 
subsonic domain but complex otherwise. The new boundary 
value problem is set by prescribing the real part of iM£,£) 
on the boundary of the unit circle. This problem has, as 
yet, not been shown to be well posed, although its numerical 
solution was found successfully. Our contribution will be 
to show that a corresponding problem for the Tricomi equation 
is well posed. 

Under the change of coordinates 



equation ( 1 . 1 ) can be written in the canonical form 
(1.2) Au + 3 ^— u y = 0 . 

We consider this equation in the four-dimensional space 
where the variables X and Y are complex. Points in the 
upper unit semicircle of the real (X,Y) -plane correspond 
to points in the elliptic half-plane of the variables x and y, 
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which we shall call the Tricomi variables. Points in the 
lower unit semicircle of the real (X,Y) -plane correspond 
to a complex surface in the space of the x and y coordinates. 

A solution u(X,Y) of (1.2), real in the upper unit 
semicircle, is continued analytically to the lower unit 
semicircle over a path which avoids the singular locus Y * 0. 
The boundary value problem in which we are interested con- 
sists in prescribing the real part of u(X,Y) on the boundary 
of the real unit circle in the (X,Y) -plane. A solution of 
this boundary value problem will be given that can be 
described formally as follows: 

Two types of solutions for equation (1.2) can be found. 
One of them is of the form 

(1.3) u x (X,Y) = Y 2/3 f(X,Y 2 ) , 

2 

where f(X,Y ) is an analytic function even in Y. The other 
is of the form 

(1.4) u 2 (X,Y) = g (X , Y 2 ) , 

2 

where g(X,Y ) is again analytic and even in Y. 

The Dirichlet problem in the upper unit semicircle with 
vanishing boundary values at Y = 0 can be solved by a solution 
of the type (1.3). The mixed Dirichlet-Neumann boundary 
value problem in the same domain, with vanishing normal 
derivative on Y = 0 and prescribed values of the solution on 
the remaining boundary, can be solved by a solution of the 
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type (1.4). By analytic continuation, these solutions can 
be extended to the lower unit semicircle. The solution u 2 
will still be real. It will be of the form 

(1.5) u 2 (X,-Y) » g(X,Y 2 ) , Y > 0 . 

But the solution u^ will be of the form 

(1.6) U 1 (X,-Y) - I Y I 2/3 e 2 ^ 3 f(X,Y 2 ) , Y > 0 

when a proper branch is chosen. 

In Chapter 5 it will be shown that the solution u of 
our boundary value problem can be expressed in a unique 
way as 

(1.7) u(X,Y) - u x (X,Y) + u 2 (X,Y) , 

where u^ and u, are of the form (1.3) and (1.4) respectively. 

A uniquely invertible relationship between the boundary values 
of u and those of u^ and u 2 is obtained from the evenness 
of f and g as functions of Y. Using (1.5) and (1.6) we 
can write this relationship as 

u (X , Y) » | Y | 2/3 f(X Y 2 ) + g (X,Y 2 ) , Y > 0 , 

^.nd 

Re{u (X, Y) } = - i | Y | 2/3 f ( X , Y 2 ) + g(X,Y 2 ) , Y < 0 . 

This shows that the boundary values of and u 2 on the 
upper unit semicircle can be chosen so that the real part 
of u satisfies the prescribed boundary conditions on the 
whole perimeter of the unit circle. 

OKU ; 

OK KOOlt Qu.Ua. / 
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AN APPLICATION 


As we have mentioned in Chapter 1, our boundary value 
problem has been motivated by the problem of calculating 
shockless transonic airfoils by the method of complex char- 
acteristics described in (9). In this chapter we give a 
description of the method, specifically of those aspects 
which motivate our problem. 

For an understanding of the problem it is important to 
find a suitable set of coordinates in which it can be stated. 
We start by expressing the equations of motion in the systems 
of coordinates that are relevant. The boundary value problem 
of reference [9] is then stated. Afterwards we give a des- 
cription of transonic paths of integration that are required. 
We consider the geometry primarily in the case of the Tricomi 
equation and the Tricomi gas. This is a hypothetical gas 
whose equation of state is such that, in pertinent coordi- 
nates, the equation for the stream function becomes the 
Tricomi equation (see (4]). 

2.1. Equations of Motion and the Hodograph Transformation 

Let us consider the equations of plane flow for an 
inviscid compressible fluid. We assume steady, isentropic 
and irrotational flow. Consider an ideal polytropic gas 
whose equation of state is given by p = p\ so the process 
is isentropic. The variables p and p are respectively the 
pressure and density, and the constant, y is the adiabatic 
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exponent of the gas. 


For steady, irrotational flow the Bernoulli law becomes 



In this equation q is the speed, c is the speed of sound, 
and c* is the critical speed, which is constant. Let u and v 
be the components of the velocity vector in the (x,y) -plane. 
The continuity equation 

(2.2) (pu) ♦ (pv) - 0 

x y 

and the irrotationality condition 

(2.3) u y - v x * 0 

permit us to introduce a stream function and a velocity 
potential ♦ . Equations (2.2) and (2.3) can then be trans- 
formed, using the Bernoulli law (2.1), into second order 
partial differential equations for or We have, for 
instance, the equation 

(2.4) (c 2 -u 2 )4> - 2uv4> + (c 2 -v 2 )<)> = 0 . 

T xx Y xy yy 

The equation (2.4) is quasilinear. It can be reduced 
to a linear equation by use of the hodograph transformation. 
The hodograph transformation consists of introducing u and v 
as the independent variables. In the hodograph plane we use 
polar coordinates defined by 

ie 

u + iv = q e 
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The eq-iation (2.4) is transformed into the linear system 


(2.5) 



M 2 -! 

pq 


9 


1 * 
p q 
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where M ■ q/c is the local Mach number. 

Equations (2.5) are known in the literature as the 

Chaplygin equations. They are of the hyperbolic type for 

2 

supersonic flow, M >1, and of the elliptic type for 

2 

subsonic flow, M < 1. The ordinary differential equation 
for the characteristics of (2.5) is 

(1-M 2 ) (dq ) 2 + q 2 (d0 ) 2 - 0 . 


Introducing characteristic coordinates £#n» we obtain the 
canonical system 


qO^ - <^ 2 -l q ► “ 0 , 

( 2 . 6 ) 

q0 n ♦ q n - 0 


in the hyperbolic domain. 

We can extend the concept of characteristics to subsonic 
flow by considering analytic extension of th^ function \p into 
the four-dimensional domain of the comp ' --x variables q and 
0, or x and y (10). The solutions of (2.6) are then the 
complex characteristics and the characteristic coordinates 
f, and n become independent complex variables. In terms of 
the complex characteristic coordinates we obtain the canonical 
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equations 


(2.7) 


’5 ' lT **t ' 


♦n ‘ iT -*n 


for $ and where 



An analytic function of a characteristic coordinate is 
a new characteristic coordinate. Thus, there is a conformal 
mapping free that can be used to transform the characteristic 
coordinates conveniently. Introducing the modified speed h, 
defined by 


( 2 . 8 ) 



we can write the equations of the characteristics as 


(2.9) 


log h + iO * n* 
log h - i0 * £. 


By use of a confoimal transformation new characteristic 
coordinates s * s(£) and t ■ t(n) can be defined as 

(2.10) s * h(q)e" i0 , t - h(q)e i0 . 

Then a further mapping s = e f ^ * , t = e f< ^ can be intro 

duced to transform the flow onto the unit circle. 

The system (2.7) can be transformed putting 

= 0 , Y = log h - j /l-M 2 ‘-~j 

c* 


( 2 . 11 ) 


X 
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into the linear equation in canonical form 

( 2 . 12 ) *xx * *w * fl? *V ■ 0 

f or , with t =* t + . After another change of coordinates 
defined by 



the Chaplygin system transforms into the equation 


(2. 12) 

K <‘”%e * 

Kc - 0 ' 

with 

2 

l-M 2 


K (0 ) * T 

■ — y- • 


This be viewed as a generalized Tricomi equation. 

"ne characteristic coordinates s and t introduced in 
(2.10) have the property that, for M < 1, s * t if and 
only if q and 0 are real. Two complex coordinates with this 
property aie called conjugate coordinates (1). More 
generally let x and y be any two complex variables, and let 

* 

z = x + iy , z = x - iy . 

We see that z* = z if and only if x and y are real. There- 
fore wj call z and z* conjugate coordinates. 

Starting with the complex coordinates x and y there s 
a simple way of generating new conjugate coordinates. Let 
f ( z ) be an analytic function which is real on the real axis. 
Define £ = f(z) and n ■ f(z*). Then £ and n are conjugate 
coordinates. For if x and y are real, then z = z , so by 
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the Schwarz reflection principle f(z*) - f(z) » f(z) 
and C ■ n. The converse is also true, provided that f is 
one to one in the domain in question. 


2.2. The Inverse Problem for Supercrit ic al Airfoil Desig n 

The use of the hodograph transformation is appropriate 
for solving inverse problems in transonic gas dynamics. By 
the inverse method a continuous flow can be computed in terms 
of hodograph coordinates. Consider the streamlines in the 
physical plane. If the streamline p - 0 is closed and not 
self intersecting, it can be taken as a profile which 
generates the flow. 

A new idea is given in [9J to solve the problem of air- 
foil design. This idea is to map the unknown domain of the 
flow in the hodograph plane conformally onto the unit circle 
of the plane of one characteristic coordinate. The speed q 
on the airfoil is prescribed as a function q(s) of the arc 
length s. What is intended is to have the equations of motion 
solved for the stream function in such a way that the boundary 
of the unit circle, for subsonic points, is a streamline and 
furthermore its image in the hodograph plane acquires the 
prescribed values of q over the profile. 

More general conjugate characteristic coordinates f, and 
n can be introduced by means of the formulas 


(2.13) 



where the map function f ( C ) is an analytic function defined 
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in the unit circle | £ \ < 1, and s and t have been defined in 
(2.10). The characteristic coordinates £ and n are conjugate 
in the subsonic domain, so points in the real subsonic domain 
correspond to points £ ■ n. Beyond the sonic line points 
£ ■ n do not, however, correspond to real supersonic flow. 

The stream function , ) is thus real in the subsonic 
domain but complex otherwise. 

A new boundary value problem was formulated in 1 9) for 
the stream function <J»(£»h). Let be a solution of the 
canonical system (2.7) that is allowed to have a pole together 
with a logarithmic singularity in the hodoqraph plane. The 
boundary condition 

Re {*(£,£)} - 0 , |£| - 1 

is imposed. For this purpose a suitable branch of the 
analytic function i(>{£,n) has to be specified. 

Paths which are complex conjugates of each other are used 
to compute the solution in the real subsonic region (1,9,14). 
Given the two initial character ist ic planes £ = £ q and 
n = n 0 , a path £ = £(o) is defined in the plane n = n 0 . 

The conjugate path n ( t ) is defined in the characteristic plane 
f. = as the set of points n = £ (o) . Consider the rectangle 
formed by taking the Cartesian product of these paths. The 
diagonal of the rectangle is given by the set of points 
(£,n) = (£,£) and they correspond to real subsonic points 
until the sonic line is reached. 
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To compute points below the sonic line special paths 

2 

have to be defined that avoid the complex sonic surface M »1, 
where the canonical equations become singular. To achieve 
this the paths have to be separated so that they are no 
longer conjugate. An analytic continuation of the solution 
around the sonic line that is appropriate for the boundary 
value problem is defined by such integration paths. 

The geometry in the space of the canonical coordinates 
X and Y, in which our boundary value problem for the Tricomi 
equation is formulated, is much simpler. In this space the 
complex characteristic coordinates are given by 


and 


C - X + iY, 
n = X - iY . 


The sonic surface corresponds to the complex plane Y = 0. 
For analytic continuation around Y = 0, a path that connects 
the points (X,Y) = (0,1) and (X,Y) = (0,-1) can simply be 
defined by 

X = 0 , Y = e Hl , 0 < 0 < tt . 


The initial characteristic planes 
(0,1) are given by C = i and n = -i. 
S-path and n-path are respectively 


iY 


tt i/2 i0 
e e 


through che point 
The corresponding 


0 < 0 < Tt , 


and 


n 


iY 


e 3ni/2 e i0 _ 


0 < 0 < IT . 
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In Figure 1 we present these paths. They make it possible 
to extend solutions of (1.2) from the upper unit semicircle 
in the X,Y plane through the complex domain into the lower 
unit semicircle. 
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3. 


GREEN'S AND NEUMANN'S FUNCTIONS AND THE POISSON FORMULAS 


In this chapter we give closed formulas for the solutions 
of the Dirichlet and mixed Dirichlet and Neumann problems 
in the upper unit semicircle for the more general equation 

(3.1) u xx + u yy + £ u y ■ 0 , 

with values of k to be specified. This is the Euler-Poisson- 
Darboux equation and is also known in the literature as the 
generalized axially symmetric potential equation [27). 

For k a positive integer (3.1) is the equation of a 
potential with axial symmetry in space of dimension n = k + 2. 
Weinstein has considered equation (3.1) as the equation 
for the axially symmetric potential in a fictitious space 
of dimension k+2, with k positive and real (25-28). The 
parameter k will be referred to as the index of the 
equat ion . 

A fundamental solution valid for values of k > 0 is 
given in Section 3.1. This fundamental solution is the 
Neumann function for the upper half-plane. In Section 3.2 
we obtain the Green's function for the upper unit semicircle 
when 0 < k < 1. In Section 3 .3 we use this Green's function 
to obtain a Poisson formula for the unit semicircle; 
that is, we give a representation of the solution with 
vanishing boundary values on Y * 0 in terms of its values 
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on the boundary of the semicircle. In Section 3.4 we obtain 
tor k > 0 the representation of the solution when its 
values are prescribed on the boundary of the unit semicircle 
and the normal derivative vanishes on the X axis. We obtain 
these two representations independently. However, a 
"correspondence principle" (2b) that will be used in 
Section 3.2 relates solutions of (3.1) with different 
indices ami makes it possible to deduce the solution of 
the Di rich let problem from that of the Neumann problem. 


) . 1 The Fundamental Solution 

A fundamental solution for (3.1) with index k > 0 is 
given by 


(3.2) S(X,Y;f.,n) 



I (X-U 


. k-1 , 

sin a da 

73 . 23 .7 r 2 a T ^ r/i 

(Y~n) + 4Yn sin j ) 


9 


where the singularity is at the point (\,n). This is a general- 
i .'.it ion of the Gauss formula for the gravitational 
potential of a circular wire in a three-dimensional space, 
to which it reduces if we make k * 1 ([10), p. 148). 

Similar formulas have been given by Weinstein (25-28), 
although there is some confusion in the literature. 
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Fundamental solutions have been known for the Tricomi equa- 
tion since the time of his original paper 123). 

It can be checked by direct calculation that (3.2) is 
a solution of equation (3.1). Moreover, as a function of 
the parameters £ and n , (3.2) is a solution of the adjoint 
equat ion 


(3.3) 


k k 

v re + v - - v + — *■ v » 0 

££ nn n n 


We introduce the notation 


(X-£ ) 2 + (Y-n) 2 , Q? - 


(X -£ ) 2 + (Y+n ) 2 , 


(3.4) 


With this notation (3.2) can be written in the form [25) 


(3.5) S (X , Y ; £ , n ) - (~) k/2 


t k-1 ( 1-t 2 ) 


(k-2)/2 r 2 ^-k/2 


[fi * 


This is a hypergeometric integral. It can be expressed as 


(3.6) S(X,Y;£,n) 


2 k ” 1 n k (r (y) ) 2 


TTFT 


, ( k k . . 

|2'2 ;k;1_l 


by use of the relationship (18) 
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1 

(3.7) F(a,b;c;z) - j t b_1 (l-t) C ' b “ 1 (l-zt)" a dt, 

0 

valid for Re c > Re b > 0 , |arg(l-z)| < tt . 

The hypcrgeometric differential equation has three singular 
points, and at each of those points one of the two linearly 
independent solutions has a logarithmic singularity. Using 
known relations for the hypergeometric function 117], we 
can rewrite formula (3.6) as 

k k 

(3.8) S(X,Y;£,n) - |F(j,|;l;q 2 ) log i + G(X,YK,n)} , 

Q 1 

where G is a regular function. The coefficient of the 
logarithmic term, namely, 


(3.9) 


A(X,Y;C,n) 


_k k 
2 n 



F 


k k . 
2 * 2 ; 1 


9 


is the Riemann function of equation (3.1), which shows 
that (3.8) has the correct singularity for a fundamental 
solution . 
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3 . 2 The Green's Function for the Unit Semicircle 

Let z(k) bo a solution of equation (3.1), where 
temporarily we let the index k bo either positive or 
negative. The following "correspondence principle" 
has been introduced by Darboux and extensively used by 
Weinstein and Diaz (7,26): 


(3.10) 


z (k) - Y 1 " k z (2-k) . 


The proof is done by direct verification. For 0 «. k < 2 
it relates solutions with positive index. Using this 
principle wo have from (3.2) a now fundamental solution (28) 


(3.11) U ( X , Y ; £ , n ) 


nY 


1-k 


1-k 

sin a da 


2 2 

l (X-U +(Y-n) +4Yn sin 


rrr 

2 


(2-JTT72 


The results of Section 3.1 applied to the solution of index 
(2-k) show that (3.11) has the proper singularity. For 
0 < k < 1 the expression (3.11) is a Green's function 
relative to the upper half-plane, since it vanishes on 
the X axis. 

We use the method of inversion to obtain the Green's 
function for the unit semicircle. Our idea is to extend 
Kelvin's transformation to equation (3.1). More precisely, 
it can bo proved by direct calculation that if u(X,Y) 
is a solution of (3.1), then so also is 
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u' (X,Y) - r" k u (X ' , Y ' ) , 

where 



Applying Kelvin's transformation to (3.11), we obtain a 
new solution of (3.1) with respect to the variables X and 
Y, namely, 

(3.12) U'(X,Y;C,n) 

11 1 -k 

_ v l-k k-2 Bin a da 

j, [ (x'-0 2 +(Y'-n) 2 +4Y'n sin 2 j j<2-k)/2* 

We recall that (3.11) is a solution of the adjoint 
equation (3.3) as a function of the variables £ and n. 

Using the fact that if u(X,Y) is a solution of (3.1), 
we find that 

(3.13) v (X, Y) = Y k u (X, Y) 

is a solution of (3.3). Applying Kelvin's transformation, 
we obtain the new solution 


(3.14) V (X, Y ; £ , n ) 


,,1-k k-2 f sin* k a da 

nY p 2 2 T~ a / 2-k ) /i • 

J 0 l(x-C') +(Y-n') +4Yn' sin^ | ] u K,/ ^ 
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of (3.1), where 


S' 




C 2 * n 2 


Formula (3.14) defines a fundamental solution of (3.1) 
with singularity at the point (C',n'), which is the 
inverse image of (£,n) with respect to the unit circle. 

It can be verified that (3.11) and (3.14) take on 
the same values on the boundary of the unit semicircle. 
For that we have just to notice that for (X,Y) on the 
circumference of radius 1, the relation 



holds, where 

R 2 - (X-S) 2 + (Y-n) 2 , R 2 * (X-C) 2 + (Y-n') 2 . 


Subtracting (3.14) from (3.11), we obtain the Green's 
function for the unit semicircle. Let us introduce the 
abridged notation 


( 3 . 15 ) I (X,Y;C,n) 

r 


l-k 

sin a da 

[ (X-') 2 +(Y-n) 2 +4Yn sin 2 £ ) p/2 


Then the Green's function for (3.1) in the unit semicircle, 
with index 0 < k < 1, is 
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(3.16) G(X,Y;4,n> - nY 1 “ k ( I ? _ R (X , Y » C , n > - . K ~ 2 l ? _ k (X , > j r ■ . • ) ! . 

In the next section we shall insert this result into 
Green's formula to obtain a Poisson formula for the unit 
semicircle. 


J . 3 The Poisson Formul a for the Unit Semicircle 

The Green's function obtained in Section 3. 2 can be used 
to solve the Dirichlet problem in the unit semicircle. 

We want to obtain the solution to (3.1) which takes given 
continuous values on the upper boundary 3D of the unit 
semicircle D and which vanishes at Y ■ 0. We start from 
the representation formula [10] 

(3.17) U(X,Y) - j! J u(4,n) G M ( f. , n ; X , Y ) ds , 

3D 

where G., is the Green's function of the adjoint equation, 
3/3v stands for the inner normal derivative , and s is the 
arc length. 

For G„ we have the interchange rule 
G M (4,n;X,Y) = G L (x,Y;4,n> , 
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where G^ in the Green's function (3.16) of (3.1). Formula 

(3.17) can be applied to any domain such that the coeffi- 
cients of equation (3.1) are regular up to the boundary. 
In particular, we can apply it to a domain D ( obtained by 
deleting the strip y < e from the unit semicircle. By 
continuity we can then extend the formula to the case in 
which D is the unit semicircle itself. We can then write 

(3.17) as 

n 


(3.18) u (X , Y ) 


1 

2" 


P k ' 2 


fU) j. | Y 1 ” k p sin ♦U 2 _ k (X,Y;C,n) 
l 2 . k (x ( Yii',n'))} d* , 


where f ($) represents the values of the solution u on the 
boundary of the semicircle, and where use has been made of 
the fact that the solution vanishes at Y * 0. In this 
context we have introduced the notation 

£ * p cos $ , n ■ p sin $ . 

The kernel in (3.18) can be computed by a straightforward 
calculation. We obtain, using the notation (3.15), 


(3.19) 


3G 

3p 


(2-k)Y 1_k sin 


<Ml 2 _ k (X,Y;S,n) 


28I 4 _ k (X,YK,n) ) 


at p = 1, where 


8 = X cos $ + Y sin $ 


2Y sin $ 


2 a 

sin 2 . 
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Manipulatinq formula (3-19), we obtain 


(3.20) ~ - (2-k) Y 1 ‘ k sin $(r 2 -l) I 4 _ R <X,Y;C,n) 


Therefore, the solution (3.18) can be written as 


(3.21) u (X , Y ) - 


2 - 


T1 

k 

v 


f («) (l-r 2 )Y 1_k 


_ ./ f sin 1 k a da 1 .. 

' k l (X-U 2 +(Y-n) 2 + 4Yn Sin 2 


s i 


when 


C * cos 0 , n » sin # . 


The hypergeometric integral under the bracket can be 
transformed into a hypergeometric function using results 
similar to those presented in Section 3.1. We obtain for 
this integral the formula 


I 4 _ R (X,Y;£,n) 


,1-k 


(rii- |)) 2 
r (2-v ) 



k 

2 ' 




)• 


Making use of the quadratic transformation 


F ( 2- 1- 2-k ; z ) = ( 1-z ) " ( 1_k/4 ] F 1- |, 


k 3 k z‘ 


4 ' 2 2 ' 4z-4 I ' 


for the hypergeometric function (118], p. 51), we obtain 
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from (3.21) after further transformations, the representa 
t ion 


(3.22) u (X , Y ) - cY 


n 0 0 

l-k f f (») [1- (X*+Y ) ) sin < 
l (1+X 2 +Y i -2X cos 


/i . * 1 k 3 k 

1 4' 2 4' 2 2' 


4Y 2 .In 2 * jj d< , 


(1+X 2 4Y 2 -2X cos <f>) 


where 


2-k 

A 


(m - j>) 
r (2-k ) 


This formula generalizes a result of Germain and Bader [8] 
that was given for the case k = ~ , and is the representa- 
tion we have been seeking. It solves the Dirichlet problem 
for the Euler- Poisson-Darboux equation with index 0 < k < 1 
in the unit semicircle and with vanishing values of the 
solution at Y = 0. 
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3.4 The Mixed Dir ichlet-Neumann Problem 


The second problem we are interested in is the mixed 
boundary value problem in which the values of the solution 
are given on the upper boundary of the semicircle and the 
normal derivative vanishes at Y * 0. We ask that the solu- 
tion be continuous up to the boundary. We have already 
found in (3.2) a Neumann's function for the upper half-plane, 
valid for k > 0. It can easily be checked that the normal 
derivative of this function vanishes at Y = 0. 

In order to obtain a representation formula for the 
solution we have to devise a formula different from (3.17), 
since we are dealing with a Neumann condition at Y * 0. 

The difficulty is that we have an equation which is not 
selt-ad joint , so we have to find a function that has simul- 
taneously the properties of a Neumann's function for our 
equation and of a Green's function for the adjoint equation. 

We start with the identity [10] 

(3.23) || (vLful -uM[v] ) dX dY + | [v ^ -u ^ +(a +b |^)uv)ds 
D 3D . 0 , 

where 

L[u] = u xx + u yy + a(X,Y)u x + b(X,Y)u y + c(X,Y) = 0 

and where M[v] =0 is the adjoint equation. 
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This formula is valid if the coefficients a, b, and c are 
regular up to the boundary of the domain D. In our case 
we have 


a * 0 , 



c * 0 . 


Formula (3.2) can be written as 

n l - n k J k (x,YK,n> , 


where we use the notation 


• k_1 a 

, iw w r Sin Cl do 

V X ' Y »' ~ 2 2 2 a p/2 

P j Q KX-C) %(Y-n) +4Yn sin^ 


Then , 


N m = Y k J R (X , Y ; C , n ) 


is a fundamental solution for the adioint equation 

Mlvl - v xx * V YY - £ v Y + tj. v - 0 

By use of the Kelvin transformat ion we obtain, as 
in Section 3.3, the Green's function 


(3.24) 


G M (X,Y;C,n) 


U k (x,Y;C,n) - 


-k 


J R (X,Y;f.' ,n’ ) 1 , 
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where 






2 .. 

P 



+ 


2 

n 


The function G M vanishes on the whole boundary of the semi- 
circle, and the expression in square brackets in (3.24), 
is a solution of L(u] = 0 which has vanishing derivative 
with respect to Y at Y = 0. 

We apply the identity (3.23) to a solution u(X,Y) of 

(3.1) and to v = G„ . We take as domain of integration D / 

m e 

the unit semicircle with a strip 0 <_ Y < e deleted. Since 
G m behaves like a Green's function on the upper boundary 
f of the semicircle, we obtain 


2iru(f;,n) = 


U(X ' V> 3^ °M dS 


-It 


9u 
M 9Y 


+ k r 

U 9Y Y UG M 


dX 


where we note by y the segment Y = e of D ( . The second 
integral gives, after a cancellation and putting 

= J R (X , Y ; £ ' ,n ' ) , 


, 3jj _ 9G M k 

J M 9Y 9Y Y UG M 


dX = 


■ 1 W 


k 3 , , -k ’ 
uy 5 y(J k -P J k ) 


dX 
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If we now let i >0, the whole inteqral over > vanishes 
once we assume .'u/i'Y bounded. Thus we obtain the formula 


u ( f, , n ) 


71 

‘ ^ I 


f(*> y r Y k (J k (X,Y;C,n) 


p" k J k (X,Y»f.' ,n’) | <U, 


or, equivalently. 


u ( X , Y ) 


1 

2rt 


f(4>) n k U k (X,Y;f.,n) 


r" k J k (x’ ,Y' K,n) ) d*. 


where 


X' * , Y' - -y- , and r 2 - X 2 + Y 2 . 

r“ 


~1 

r 


Followinq the same steps as in Section 3.3 we are led 
to the formula 


u (X , Y ) 


II 

k 

2? 


. k . 

sin ^ 


f (*) 


2, . k-1 . ^ 

(1-r ) sin a da 1 . . 

777TT2I73 — r27^ — •' ' 2 a. ~nr+T r/?/ cS4 ' ' 

l(X-f.) +(Y-n) +4Yn sin ,] 


where 


f, ,* cos $ , n * sin $ . 

This formula can be written, by usinq similar relations for 
the hype rgeome trie function, as 


2H 



(3.25) u (X , Y ) 


c 


. 


f(*) ( 1“ (X + Y 2 ) ) 

2 — 2 

(1+X +Y -2X cos *) 


• k. 

sin * 

TkT2T7? 





4Y 2 sin 2 * 

5 — 5 5 ~J d, P' 

(1+X +Y -2X cos *) 


with 

k-2 (n£)) 2 
c •- 


The result is valid for k > 0. It is the desired represen- 
tation of the solution of the mixed Dirichlet-Neumann problem 
in the unit semicircle, with vanishing normal derivative on 
the X axis, in terms of the values of the solution on the 
boundary of the semicircle. 

The correspondence principle (3.10) can be applied to 

the solution (3.25) with index k' = 2-k and boundary values 
k- 1 

sin * f(*) , 0 < k < 1. We obtain from it the solution 

(3.22) of the Dirichlet problem with boundary values f(*), 
as was mentioned earlier. 
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4. 


THE REFLECTION PRINCIPLE 


We establish in this chapter reflection laws for 
solutions of the general Euler-Poisson-Darboux equation 

(4.1) u xx u yy + £ u y « 0 , 0 < k < 1 , 

which make it possible to extend them from the upper to the 
lower unit semicircle. They will be used in Chapter 5 to 
prove our main theorem about the transonic boundary value 
problem. However, they have an interest of their own. 

General theorems about reflection laws for solutions of 
elliptic partial differential equations with analytic condi- 
tions on the boundary are well known (Garabedian (11), Lewy 
(16)). The special feature of our case is that the boundary 
over which the reflection is going to be made is a singular 
line for the equation. 

It is known (Hyman (12)) that a solution of (4.1) 
which is analytic in some neighborhood of a segment of the 
X axis is necessarily even in Y and is uniquely determined 
by its values on this segment. This is a local theorem 
based on power series expansion of the analytic solution, 
valid for k different from 0 or a negative integer. The 
result for k = 1 is classic (13). However, this theorem does 
not apply to solutions with a branch singularity of the type 


30 


origins. 
qp poor, qualu 


1 


(4.2) u(X,Y) - Y 1_k f(X,Y 2 ) . 

Our reprosentat ion formula (3.22) is an example of this type 
of solution. We shall derive reflection laws that are valid 
in the large and permit one to reflect solutions of the type 

(4.2) for which the X axis is a branch line. We use analy- 
tic continuation of the solutions into the complex domain 

of the variables X and Y to do this. 

Let the variables X and Y be extended to complex values 

X = X x + ix 2 , Y « Y l + iY 2 = |Y|e 11 ' . 

We define the path > by 

Y ■ j(X,Y: X * 0, Y = re 1 ', 0 <0 < it 

where 0 - r < 1. Consider the analytic continuation of the 
solutions (3.22) and (3.25) along this path y* The explicit 
form of the solutions guarantees that the continuation can 
be performed. Two points A and B, with A in the upper unit 
semicircle and B in the lower, can bo connected by a path of 
the form 

(4.3) T = Yj u Y u Y 2 , 

where Yj is a path contained in the upper semicircle which 
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connects the points A and (0,r), and *- 8 the corresponding 
path which connects B and (0,-r) in the lower semicircle. 

By the monodromy theorem, any other path which does not wind 
around the X axis will give the same result. 

We obtain in this fashion the following theorem: If 
u(X,Y) is a solution to the Dirichlet problem for equation 
(4.1) in the upper unit semicircle with zero boundary values 
u = 0 at Y ■ 0, then it can be extended to the lower semicircle 
by the reflection law, 

(4.4) u (X , -Y ) = e (1_k) 1 u (X ,Y) , Y > 0 , 

where analytic continuation has been made along the path 
(4.3). Similarly, the solution of the mixed Dirichlet and 
Neumann problem in the same domain, with vanishing normal 
derivative u y = 0 at Y = 0, can be extended to the lower 
semicircle according to the law 

(4.5) u(X,-Y) = u (X,Y) , Y > 0 . 

Clearly, the reflection rules depend only on the form of the 
boundary condition at Y = 0 and are valid in more general 
domains . 
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5. SOLUTION OF THE BOUNDARY VALUE PROBLEM 

We can now state and prove our main theorem. Consider 
the Tricomi equation 


(5.1) 


yu + 
1 xx 


u 


yy 


o , 


and the transformed equation 


(5. 2) 


u 


XX 


U YY + 


1 

3Y 


u 


0 , 


in canonical coordinates X and Y. We introduce the class of 
analytic solutions of equation (5.1) which are single-valued 
in a neighborhood of the line y = 0. Let u(X,Y) be a 
real solution of (5.2) in the upper unit semicircle 


D + ■ { 


(X,Y) : X + Y < 1 , Y > 0 


belonging to this class. Then in the complex domain of the 
extended variables X and Y the solution u(X,Y) can be continued 
analytically along a path, as described in (4.3), into the 
lower unit semicircle, 

D = ( (X,Y) : X 2 + Y 2 < 1 , Y < 0 



We shall prove that the boundary value problem 
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Re | u (X, Y) | « F ( 0 ) , X - cos 0 , Y - sin 0 , 0 < 0 < 2n , 

where F(0) is periodic and continuous, is well posed within 
the specified class of solutions u(X,Y). 

We give first an explanation of the choice of the class 
of functions that we have made. If u(x,y) is single valued 
and analytic, it admits a double power series expansion in a 
neighborhood of the origin of the form 


u(x,y) = l 


n m 

a x y 
nm 1 


Therefore, we can write 


u(x,y) * 
or 

u (X, Y) * v : 


3 ~ 3 2 - 3 

v^ (x ,y ) + yv 2 (x,y ) + y v 3 (x,y ) , 

(X , Y 2 ) + Y 2/3 v 9 (X,Y 2 ) + Y 4/3 v 3 (X,Y 2 ) 


# 


where v^ , v 2 and v 3 are regular. 

If we let the variable Y wind around the X-axis, through 
an angle of 2v or 4n, we obtain new solutions 


u' (X , Y ) 


v L (X,Y ) + 


e 4rx/3 y 2/3 V2 ( X ,y 2 ) + e 8 " l/3 y 4/3 v 3 (X , Y 2 ) 


and 


u " ( X , Y ) = v x (X,Y 2 ) + e 8TTl/3 Y 2/3 v 2 (X,Y 2 ) + e 1 6n l/3 Y 4/3 v 3 (X , Y 2 ) 


of equation (5.2), as the equation itself remains unchanged. 
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We apply now the differential operator L defined by equation 
(5.2) to u, u' and u". This gives 


L(v 1 l 

L(v 1 J 

LlVj^) 


+ L|y 2/3 v 2 l + liv 4/3 v 3 i - o 


♦ e 4 ’ 1 ' 3 LIY 2 / 3 V2 , + c 8,i/3 l[v 4/3 V3 , . „ 


♦ e 8"i/ 3 l[v 2/3 V2 , + e 16-i/3 l[v 4/3 V3) . 0 


Because the determinant of this linear system of equations 

is different from zero, it follows that each of the terms 
3 3 2 3 

v x (x,y ), yv 2 (x,y ), and y v 3 (x,y ) is a solution of the 
Tricomi equation separately. By the uniqueness theorem for 
the Cauchy problem for (5.1), we can conclude that 


2 

y v. 


(x,y 3 ) 


0 . 


Therefore, u (X,Y) admits a unique decomposition of the form 
(5.3) u (X , Y ) = V;l (X,Y 2 ) + Y 2/3 v 2 (X,Y 2 ) . 


We proceed now to prove the main theorem. Let the 
boundary function F(0) be defined by 


u (X , Y ) = h (X) ? Y > 0 , 


X 2 + Y 2 


1 
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and 

Re ju (X , Y) | - k (X) ; Y < 0 , X 2 ♦ Y 2 « 1 . 

Let u^(X,Y) be a solution of the Dirichlet problem in the 
upper unit semicircle such that 

Uj (X , Y ) » f(X) , X 2 + Y 2 - 1 , Y > 0 , 

Uj (X , 0 ) - 0 , 

with f (X) continuous and vanishing at X ■ + 1. Let u^(X,Y) 
be a solution of the mixed Dirichlet-Neumann problem, in 
the same domain, with 

u 2 (X,Y) = g (x) , 

3u_ 

(X»0) * 0 • 

Both u^lXjY) and u 2 <X,Y) can be computed explicitly by means 
of the representation formulas (3.22) and (3.25) of Chapter 3. 
By the reflection laws obtained in Chapter 4, we have 

UjfX^Y) = e 2ll/3 f(X) , Y > 0 , X 2 + Y 2 = 1 , 

and 

u 2 (X,-Y) = g (X ) , Y > 0 , X^ + Y 2 * 1 . 
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If we choose f and g so that 


f (X) * g (X) - h (X) 
f (X ) cos j ♦ g(X) ■ k(X) , 

that is, if we set 

f (X ) - | (h (X) - k (X) ) 

and 

g (X) - 5 k (X) , 

then 

(5.4) u (X , Y ) - u^X.Y) + u 2 (X,Y) 

is a solution of (5.."J which satisfies the prescribed 
boundary conditions. 

The uniqueness of the solution of the boundary value 
problem follows immediately. The representation formulas 
(3.22) and (3.25) show that Uj (X,Y) and u^(X,Y) have the 
form 

u x { X , Y ) = Y 2 / 3 Vj (X,Y 2 ) 

and 

u 2 (X,Y) - v 2 (X,Y 2 ) 

respectively. The decomposition (5.4) is therefore of the 
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form (5.3) and is unique. The boundary condition on u is 
seen to impose Dirichlet and Dirichlet-Neumann data on u^ 
and u^ respectively, which determine them uniquely, complet- 
ing the proof. The continuous dependence on the boundary 
values is also a consequence of the explicit representation 
we have given of the solution. 

We remark that the theorem could be extended to arbi- 
trary domains symmetric with respect to the X axis by 
performing reflection in the manner described in this work, 
but resorting to a more general treatment of the Dirichlet 
and Dirichlet-Neumann problem for Uj and u 2 . 
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